Theory of exchange coupling in disordered magnetic multilayers 
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■ We consider mechamism of exchange coupling based on interaction between 

^ ; 

^ ■ electrons in nonmagnetic layer. Depending on ratio of inverse time of diffu- 

I; 

^ ■ conducting electrons this mechanism describes transition from ferromagnetic 

Ch ■ 

Q ■ to concollinear ordering of magnetizations of ferromagnetic layers. 

^ : I. INTRODUCTION AND MAIN RESULTS. 

^ ■ 

00 ; 

^ I In metallic ferromagnet-nonferromagnet-ferromagnet multilayers (see fig. 1) magnetic 

o : 

CN I structure oscillates between ferromagnetic and antiferromagnetic orientations of the ferro- 
O . 

+^ I magnets' magnetizations as a function of thickness of nonmagnetic metal L with a period of 

S I order of the Fermi wave length [HI [§] @] 0- The explanation of this phenomenon is 

' based on the fact that the interlayer coupling is due to Ruderman-Kittel interaction between 

o : 

H ' electron spins in different ferromagnets. 

Further investigations discovered structures with perpendicular orientations of the fer- 



I romagnets' magnetizations (see for review 0). Often phenomenological coupling between 
magnetizations of ferromagnetic layers can be represented as sum of bilinear and biquadratic 
contributions 

E {if) = Ji cos + J2 COS^ Lf (1) 
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Here ^ is angle between directions of magnetizations of ferromagnetic films. Bilinear 
constant Ji oscillates as function of interlayer distance L. In case of large positive biquadratic 
constant J2 minimum of E {ip) corresponds to = 7r/2 . As explained by Slonzevskii large 
positive biquadratic coupling might be result of spatial fluctuations of bilinear coupling Ji 
due to ferromagnet-nonferromagnet surface roughness [|^. 

In disordered system, when L is larger than electron mean free path / , RKKY interaction 
(Ji), averaged over realizations of scattering potential exponentially decreases [If. At the 
same time fluctuations of local exchange become much larger than (Ji) giving rise to 
biquadratic contribution J2 >> \ {Ji)\ HlQj . 

Here we propose mechanism of coupling in disordered multilayers based on interaction 
between electrons in nonmagnetic layer. Spin fluctuations in system of interacting electrons 
give rise to contribution to thermodynamic potential which depends on magnetic field 



or, in our case, on relative orientation of magnetizations in ferromagnetic layers. Here we 
show that in magnetic multilayer this mechanism describes transition between ferromagnetic 
and noncolinear ordering with increasing distance between ferromagnetic layers or value of 
ferromagnetic splitting of conducting electrons. 

We assume that magnetic multilayer can be described by Hamiltonian 

H = Ho + ee^c J dr^^ (r) n (z) aap^(3 (r) + H^nt (2) 

F 

Here Hq is Hamiltonian of free electrons in random field. Second term describes exchange 
field in ferromagnetic layers, eexc is ferromagnetic splitting of conducting electrons, n (z) is 
unit vector of direction of magnetization of ferromagnetic layers. n{z) = rii at 2; < — L/2 
and n (z) = n2 at z > L/2 as it is shown on figure 1. (r) and (r) are creation 
and annihilation operators, cTq,^ are Pauli matrixes. Integration in second term is over 
ferromagnetic layers. The last term Hint describes Coulomb interaction between electrons in 
nonmagnetic layer. We assume that interaction in ferromagnetic layers is taken into account 
self consistently in e^xc- 

Details of calculation are given in the last part of the paper. Here we present the main 



results. Characteristic energies in the problem are ferromagnetic splitting of conducting 
electrons e^xc and Thouless energy = D / . D is diffusion constant of conduction 
electrons. We assume that it is the same in nonmagnetic and ferromagnetic layers. 

In case of small thickness when Ec > eexc coupling between ferromagnetic layers has 
bilinear form and coupling energy per unit area is 



Here F is characteristic constant of interaction in diffusion channel [|rT|. It is positive for 



Coulomb repulsion between electrons. Let us note that in this regime coupling (3) does not 
depend on L . Minimum of (3) corresponds to ferromagnetic orientation of magnetizations 
in multilayer ip = 0. Note that result is obtained in limit when L > I , or e^xc smaller that 
inverse mean free time D/P. 

At larger distance L when Ec < eexc coupling has biquadratic form and coupling energy 
per unit area is 

p 

^(<^) - -^^^^^cCosV (4) 

This quantity decreases as with increasing distance. Minimum of coupling energy 
corresponds to noncollinear state ip = ti/2. 

Both expressions are given for the case of infinite thickness of ferromagnetic layers. 



Calculation show that in case d> J D / eexc coupling weakly depends on d. 



Results (3) and (4) are also valid provided L < ^ D /T . In opposite case coupling energy 

— W ^ L 

Let us compare results (2) and (3) with biquadratic contribution due to mesoscopic 
fluctuations of RKKY interaction which is J2 ~ 'J2~^ a-t Ee < eexc and t/j ~ AdE 



Ec > eexc- Here A is an intralayer ferromagnetic stiffness and thickness d > yD/eexc ■ 

The quantity J2 decreases with L much faster than (4). Also for eexc/ Ad << 1 , -F ~ 1 
coupling energy given by expressions (3) and (4) is larger than biquadratic contribution due 
to mesoscopic fluctuations of RKKY in whole range of distances. In this case with increasing 



distance L » I system is undergo transition between ferromagnetic and perpendicular 
if = Tx/2 ordering. Such transition was observed in |T2|. 



II. DERIVATION OF RESULTS. 



Correction to thermodynamic potential which depends on eexc^ {z) is given by expression 
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Here constant F describes screened Coulomb interaction in diffusion channel. u;„ = 2TmT 
is Matsubara frequency, r is electron mean free time. 
Diffusion ladder satisfies equation 



= 5{Z- Z') SafsS^rj 

It is convenient to present solution of equation (6) at 1^1 < L/2 in the form 



(6) 



= i-Q') + exp (Qz) + '-^^^-^^-^^c.,S,, (7) 



Here we introduce Q = y q"^ + U is matrix of relative rotation of magnetizations 

of ferromagnetic layers. In case when direction of magnetization in ferromagnetic layer 
z < —L/2 is directed along z axes n (2;) = (0,0,1) and at z > L/2 direction is n (2;) = 
(siny9, 0, cosy^), it is matrix of rotation along y axes U = exp (^~^^y)' 

For simplicity we consider limit of semiinfinite ferromagnetic layers. More detail con- 



sideration shows that at d > yD/e^xc results weakly depend on thickness of ferromagnetic 
layers. It is convenient to introduce boundary conditions for diffusion ladder at ferromagnet- 
nonferromagnet interfaces taking into account that according to equation (6) in coordinate 
system where spins are directed along magnetization, components of ladder with antiparallel 
spins decreases as exp {—Qi \z\) and exp {—Ql \z\) at \z\ > L/2, where Qi = J q"^ + ^ ~^ ^^^'^ 
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. Components of ladder with parallel spins decreases as e:K.p{—Q\z\) at \z\ > L/2. At 
z — —L/2 where n {z) = (0, 0, 1) boundary conditions are 



d 
dz 



Here we introduce projectors of spins on 2;-axes P± — — — ^. 

The same kind of boundary conditions can be introduced for rotated diffusion ladder 
Ua-yD/j^U^ at z — L/2 . Solving system of equations (6,8) we obtain 



\ui„t\<1 ^ 



X 



[(|Ap-(ReA)2)cosv5+(lA|''-i(RcA)2)cos2.^-i(RcA)2+RcA(l-i|A|2-ilApcos2 <p) 



(9) 



[l-(ReA)2+2(|A|2-(ReA)2)cos¥'+(|A|'*-(ReA)2)cos2 



Here A = — ^ exp {—QL). Expression (9) contains divergent terms, which do not 

[Qi + Q) 

depend (p and must be subtracted. 

In limit of large exchange splitting when \Qi\ > Q parameter A = exp(— Qi^) is real. 
In this case energy is function of cos^ (p. Subtracting in expression (9) terms which do no 
depend on angle we obtain 

F ^ , , r d\ L(l-A^) A^eosV 



Main contribution in expression (10) is from region where A < 1 , denominator therefore 
gives only small correction. Neglecting it we obtain expression (4). 

In opposite limit of small exchange splitting Re A ~ , |A| ^ and to the order e^^ 
coupling energy is proportional to cos (p. 

5n{^)^-^T E l^"l/7^^^°«^ (11) 

^ \^nT\<l;a,(3 (27r) L)Q 

Calculating (11) at zero temperature we obtain (3). Transition between limits (10) and (11) 
occurs at egx ~ D/L'^. 
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